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Abstract 

In  this  paper,  operating  benefits  from  demand- side  load  management  are  evaluated  for  a  Proton  Exchange  Membrane  (PEM)  Fuel  Cell  Power 
Plant  (FCPP).  For  reliability  modeling  and  evaluation  of  the  PEM  FCPP,  a  state- space  generation  model  for  a  stand-alone  PEM  fuel  cell  that 
calculates  the  system  availability  and  the  expected  energy  not  supplied  (EENS)  index  has  been  developed.  A  systematic  technique  and  detailed 
computer  simulation  software  for  a  stand-alone  PEM  fuel  cell  station  reliability  assessment  have  been  built.  The  suggested  technique  can  be  used 
for  practical  engineering  applications  to  provide  information  for  stand-alone  FC  generating  station  planning,  design,  and  operation.  The  simulation 
results  are  obtained  using  the  MATLAB  software  for  a  5  kW  stand-alone  PEM  fuel  cell  that  supplies  a  typical  residential  house. 

©  2005  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

Considerable  attention  has  been  devoted  to  distributed 
sources  of  energy  for  meeting  the  power  demand  instead  of 
constructing  new  conventional  power  plants  due  to  better  power 
quality,  reliability,  portability  and  ecological  constraints.  Among 
the  various  types  of  distributed  generation,  fuel  cells  (FCs),  par¬ 
ticularly  PEM  FCPPs  generated  tremendous  interest  for  electric¬ 
ity  and  heat  generation  due  to  their  low  operating  temperature, 
fast  start  up  characteristics,  and  ecological  constraints  [1].  The 
use  of  Fuel  Cell  Power  Plants  (FCPPs)  is  expected  to  become 
more  widespread  in  the  near  future,  in  spite  of  their  high  current 
capital  cost.  During  normal  operation,  FCPPs  are  subjected  to 
a  number  of  possible  outage  and  derated  modes  due  to  partial 
or  full  failure  of  system  components  as  well  as  degradation  of 
the  fuel  cell  stack  and  battery.  Some  FC  failure  modes  include 
membranes  drying  out  (insufficient  humidification),  overheat¬ 
ing,  passages  clogging  up  with  water  and  freezing  of  water  in 
humidification  channels.  For  these  reasons,  fuel  cells,  particu¬ 
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larly  grid-independent  FCPPs  may  encounter  lack  of  generation 
in  meeting  the  peak  demand.  One  way  to  overcome  this  prob¬ 
lem  is  to  use  load  management  system  in  stand-alone  FCPPs. 
Load  management  in  this  study  refers  to  any  load  curtailment 
activity  that  shaves  the  peak  so  that  the  power  demand  stays  in 
the  range  of  FCPP  supply  limits.  This  action  is  also  needed  to 
increase  competitiveness  and  market  value  of  the  FCPP.  To  show 
the  impact  of  load  management  on  FC  operation,  it  is  important 
to  analyze  FCPP  system  reliability  for  likely  operational  and 
loading  conditions.  However,  very  little  information  is  avail¬ 
able  on  the  reliability  of  FCPPs  possibly  due  to  the  early  stage 
of  FCPP  technology  dependent,  insufficient  data  and  uncertain¬ 
ties.  Consequently,  the  results  presented  in  this  paper  provide 
key  information  for  individual  FCPP  planning,  design  and  oper¬ 
ation.  The  reliability  studies  related  to  FCPPs  are  generally 
investigated  from  the  viewpoint  of  improving  the  continuity 
of  electricity  supply  rather  than  the  reliability  level  of  FCPP 
itself.  Smith  and  Giancaterino  [2]  utilized  the  PEM  FCPP  in  a 
telecommunication  back-up  power  system  to  maintain  high  sys¬ 
tem  reliability.  Another  study  by  [3]  considers  the  fuel  cells  as 
a  distributed  power  generation  with  high  reliability,  durability 
and  environmental  benefits.  In  this  study  [3],  a  component-based 
state-space  representation  technique  is  considered  as  a  basic 
model  in  the  stand-alone  PEM  FCPP  reliability  evaluation.  The 
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state- space  method  includes  different  states  such  as  operation, 
derated,  fully  faulted  or  maintenance,  since  FCPPs  are  subject 
to  a  number  of  possible  outage  and  derated  states.  The  aging 
failures  in  the  state- space  generation  model  of  the  PEM  FCPP 
are  estimated  based  on  the  assumption  that  the  failure  and  repair 
rates  of  the  components  change  with  the  operation  age  of  FCPPs. 
The  functional  relationship  between  the  age  of  FCPP  and  tran¬ 
sition  rates,  namely  failure  and  repair  rates  is  estimated  based 
on  the  fuzzy  set  theory  and  expert  knowledge. 

2.  Background  and  notations 

In  a  fuel  cell,  hydrogen  is  fed  at  the  anode,  oxygen  is  fed  at 
the  cathode,  and  an  electrolyte  is  sandwiched  between  the  two 
electrodes  for  conveying  ion  e~  from  the  anode  to  the  cathode. 
Electrons  are  carried  to  the  cathode  through  both  anode  and  a 
conducting  wire,  and  a  load  is  placed  in  between.  There  are  many 
auxiliary  devices  needed  to  operate  the  FCPP,  which  take  part  in 
the  gas  and  electricity  management  and  are  used  for  regulating 
the  parameters  such  as  reactant  flow  rate,  total  pressure,  reac¬ 
tant  partial  pressure,  temperature  and  membrane  humidity  at  a 
desired  value  to  ensure  that  FCPPs  can  run  smoothly  without 
getting  the  stack  either  flooded  or  drying  out  [4] .  Accordingly, 
any  malfunctioning,  performance  loss,  and/or  failure  in  these 
auxiliaries  can  reduce  the  overall  performance  of  the  FCPP. 
This  section  introduces  an  overview  and  the  terminology  used 
throughout  this  paper. 

2.7.  PEM  fuel  cell  model 

Fuel  cells  basically  convert  chemical  energy  of  hydrocar¬ 
bon  fuels  directly  into  dc  form  of  electrical  energy.  A  FCPP 
mainly  consists  of  a  fuel-processing  unit  (reformer),  FC  stack 
and  power-conditioning  unit.  The  FC  uses  hydrogen  as  input 
and  produces  dc  power  at  the  output  of  the  stack.  A  simple  rep¬ 
resentation  of  a  FC  system  is  shown  in  Fig.  1. 

The  performance  of  a  FC  is  generally  characterized  by  using 
the  polarization  curve,  which  is  a  plot  of  the  FC  output  voltage  as 
a  function  of  load  current.  The  polarization  curve  is  computed  by 
using  the  Tafel  equation  [5],  which  subtracts  the  various  voltage 
losses  from  the  open  circuit  dc  voltage,  and  is  expressed  as 


Fstack  —  Vopen  Tohmic  ^activation  ^concentration 


Fig.  1.  Basic  fuel  cell  components. 


where, 
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In  the  above  equations,  No  is  the  cell  number,  Vo  the  open  cell 
voltage,  R  the  universal  gas  constant,  T  the  temperature  of  the 
fuel  cell  stack,  Fthe  Faraday’s  constant,  Pu2  the  hydrogen  partial 
pressure,  Ph2o  the  water  partial  pressure,  Po2  the  oxygen  partial 
pressure,  Po  the  standard  pressure,  a  represents  the  charge  trans¬ 
fer  coefficient  of  the  electrodes,  I^c  the  current  of  the  FC  stack, 
Aim  the  limiting  current  of  FC  stack,  7o  the  exchange  current 
of  FC  stack  and  c  is  the  empirical  coefficient  for  concentration 
voltage.  The  steady  state  voltage  for  one  cell  (No  =  1)  and  power 
versus  cell  current  density  is  obtained  based  on  Eq.  (1)  and  as 
shown  in  Fig.  2. 

Various  auxiliary  components  such  as  air  compressors, 
pumps,  humidification  equipment,  blower  and  coolers  are  used 
in  the  FCPP  that  are  related  to  thermodynamics  and  flow  control. 
Besides  the  power-conditioning  unit  (dc/dc  converter  plus  dc/ac 
inverter),  control  electronics,  energy  storage  and  transformers 
are  used  in  power  conversion  and  overall  system  control.  Fig.  3 
shows  a  PEM  FCPP  block  diagram  that  shows  the  auxiliary 
components  along  with  input  and  output  signals. 


2.2.  State-space  based  generation  unit  reliability 

Most  generation  units  require  auxiliary  equipment  and  there¬ 
fore  they  are  subjected  to  different  possible  derated  capacity 
states  based  on  factors  such  as  outage  of  auxiliaries,  fuel  quality 
and  environmental  conditions  [6].  Hence,  the  state-space  model 
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Fig.  3.  PEM  fuel  cell  system  block  diagram. 


of  the  generating  unit  may  include  different  states  such  as  up, 
derated  and  completely  down.  The  general  state-space  model  of 
a  generating  unit  for  the  aforementioned  case  is  given  in  Fig.  4, 
where  X  is  the  failure  rate  and  fi  is  the  repair  rate. 

In  Fig.  4,  the  indices  i=  1,2,...,/:  represent  the  number  of 
components  that  cause  unit  failure;  the  indices  7=1,2,  . . .,  m 


represent  the  number  of  components  that  affect  derating  level 
of  the  unit.  The  state-space  based  reliability  calculation  is  per¬ 
formed  using  Markov  models.  The  state-space  equation  of  Fig.  4 
can  be  written  as 


d 
d  t 


P(t )  =  A  •  P(t) 


DOWN  STATES 


DERATED  STATES 


I 


Fig.  4.  The  state-space  model  of  the  generating  unit. 
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where  A  is  the  transition  matrix  defined  as 


n 
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7-rci 

i=2 

7-12 

/  i 

—  7.2 i 

Lnl 

i=l,ifz2 

7-1  n 

'Lin 

n  —  1 

i=\ 

and 

Pit)  =  [Pi(t)  Plit)  Pi(t)  PA(t)f.  (8) 

In  Eqs.  (7)  and  (8),  Xij  represents  the  transition  rate  between  the 
ith  state  and  the  jth  state,  and  Pft)  represents  the  probability  of 
the  ith  state. 

2.3.  Load  management  system 

The  primary  goal  of  the  load  management  system  is  to 
reduce  the  overall  power  consumption  while  leveling  the  load 
by  smoothing  the  peaks  and  valleys  in  the  power  demand 
profile  of  the  system.  The  load  management  and  control  sys¬ 
tem  includes  data  acquisition  system,  user  interface,  and  smart 
energy  management  control  (SEMaC)  software.  The  SEMaC 
software  makes  decisions  regarding  the  control  actions  needed 
to  keep  the  maximum  load  below  the  FCPP  rated  capacity  at  any 
given  instant  as  shown  in  Fig.  5. 

Power  management  decisions  must  be  made  in  the  context 
of  minimizing  user  inconvenience  and  maximizing  user  comfort 
while  reducing  the  overall  power  consumption.  An  effective  sys¬ 
tem  is  virtually  unnoticeable  to  the  user,  but  results  in  substantial 
savings  in  the  annual  power  bill  by  shaving  peak  demands  in  the 
system  load  profile.  Peak  loads  are  managed  based  on  the  load 
rescheduling  algorithms.  The  load  management  system  scans 
the  data  acquisition  system  and  calculates  the  power  usage  of 
all  loads.  If  the  power  is  over  a  predefined  level,  the  peak  load  is 
managed  by  sending  an  appropriate  signal  to  the  SeMAC  system 
to  reduce  or  turn  off  the  load. 


Fig.  5.  Load  profile  comparison  with  and  without  load  management  system. 


3.  PEM  FCPP  reliability  model 

The  FCPP  includes  fuel  cell  stack  and  auxiliary  equipments 
such  as  air  compressors,  pumps,  humidifiers,  coolers  and  control 
electronics,  as  shown  in  Fig.  3.  During  normal  operation,  some 
FCPP  components  such  as  compressor,  fans,  pumps,  motors, 
temperature  and  humidity  sensors,  relays  and  other  control  elec¬ 
tronics  may  contribute  to  system  failure  or  derated  mode  for 
different  reasons  such  as  ignition  of  any  leaking  hydrogen, 
material  fatigue,  wear  outs,  break-downs,  membrane  drying  out, 
overheating  and  freezing  of  water  in  channels  [7] .  If  any  of  these 
components  go  beyond  the  operating  limits,  the  output  of  that 
component  will  be  reduced  by  a  factor.  For  instance,  insufficient 
circulating  coolant  flow  due  to  the  coolant  water  pump  failure 
may  cause  a  reduction  in  the  nominal  output.  Other  failures  may 
cause  either  a  reduction  in  nominal  output  power  or  total  sys¬ 
tem  outage.  The  results  of  the  aforementioned  auxiliary  failures 
are  considered  herein  rather  than  the  failure  modes  themselves 
to  calculate  the  effects  of  sub- system  performance  reduction  on 
overall  FCPP  performance. 

3.1.  Cooling  system 

If  the  cooling  is  insufficient  due  to  loss  of  performance  in 
the  coolant  system,  the  operating  temperature  of  the  FCPP  will 
increase  to  T  =  T  +  AT,  where  T  denotes  the  new  operating 
temperature.  An  increase  in  temperature  leads  to  reduction  in 
FCPP  output  voltage  due  to  the  fact  that  Gibbs  free  energy 
changes  inversely  with  temperature  as  shown  in  the  following 
equation. 

Pstack  —  Pstack  ~  AV 

The  other  effects  of  temperature  change  include 

•  Resistive  voltage  loss,  Pohmic  tends  to  increase  at  higher  tem¬ 
perature  due  to  dry  out  mode  of  the  FC  stack. 

•  Activation  voltage  loss,  Pactivation  increases  due  to  the  constant 

DT 

of  in  the  equation. 

•  The  only  reduction  takes  place  at  concentration  voltage  loss 

DT 

due  to  the  negative  coefficient  in  PCOncentration  expression. 
However,  percentage  voltage  loss  reduction  in  PCOncentration 
will  be  smaller  than  a  percent  voltage  loss  increase  in  Pactivation 
due  to  coefficient  a  that  takes  a  value  between  0  and  1  (gen¬ 
erally,  a  =  0.5).  Hence,  a  decrease  in  the  maximum  power 
level  (P  =  P  —  AP)  will  be  observed  for  higher  tempera¬ 
ture,  where  the  power  change  A P  is  expressed  as  A P  =  AP 

Clc  • 

Consequently,  if  the  temperature  increases  due  to  insufficient 
cooling,  the  power  supplying  capacity  of  the  FCPP  also  reduces 
as  shown  by  the  following  expression 

P  -  AP  P 

RPC[%]  = - - -  x  100  =  -  x  100  (7) 

where  RPc  is  the  percent  reduction  in  power  supplying  capacity 
of  FCPP  due  to  insufficient  cooling.  If  the  FCPP  cooling  system 
completely  fails  or  if  such  failure  leads  to  instability  in  FCPP 
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operation,  then  an  emergency  stop  function  is  activated  to  shut 
off  the  system,  and  the  system  will  go  to  off  state. 

3.2.  Humidification  system 


usage  can  be  calculated  as  follows: 


<?h2  ~  A<7h2 

m2 


x  100  = 


<7h2 


x  100 


(10) 


PEM  fuel  cells  are  quite  sensitive  to  the  humidity  of  reac¬ 
tant  gases.  Any  excess  liquid  water  from  the  humidification 
process  that  flows  into  the  FCPP  can  reduce  the  output  power 
capacity  and  may  cause  stack  damage.  The  system  performance 
loss  is  mainly  caused  by  the  internal  ohmic  resistance  increase 
resulted  from  drying-out  of  the  membrane  (insufficient  humidifi¬ 
cation).  After  the  membrane  dry-out  mode,  the  magnitude  of  the 
FCPP  stack  impedance  increases.  This  change  can  be  expressed 
by  R^c  =  Pfc  +  AT?.  However,  the  reduction  in  FCPP  output 
power  depends  on  the  severity  of  insufficient  humidification.  The 
resistive  voltage  loss,  V0hmic  increases  in  proportion  to  FCPP 
stack  resistance,  which  leads  to  output  voltage  reduction  and 
consequent  power  reduction. 

Since  humid  reactants  are  a  mixture  of  dry  gasses  and  vapor, 
dry  partial  pressure  is  the  difference  between  total  pressure  and 
the  vapor  pressure,  which  can  be  expressed  as  pj  =  pj  —  pw. 
Since  humidity  ratio  (0)  is  expressed  as  the  ratio  of  partial  pres¬ 
sure  of  water  vapor  (py)  to  saturation  pressure  (ps),  then  the 
vapor  pressure  (pv)  can  be  written  as  py  =  (p-ps.  If  Nernst  voltage 
(V0pen)  is  rewritten  in  terms  of  the  total  system  pressure  (0>t), 
we  get 


Vopen  —  E  T 


RT 

2F 


+ 


RT 
4 ~F 


In  (pr) 


where  w,  x,  and  y  are  the  constants  related  to  molar  masses 
and  concentrations  of  H2,  O2  and  partial  pressures,  denoted  as 
pn2  =  w  -  pt ,  po2  =  x  •  pr,  pu2o  =  y  •  Pt-  From  Eq.  (8),  it  is 
evident  that  if  the  total  system  pressure,  pj,  reduces  to  pj,  then 
the  open  circuit  voltage  and  the  corresponding  output  voltage 

r>  rT' 

will  be  reduced  by  ^  ln(/?x  —  pr)- 


3.3.  Fuel  supply  system 

For  the  PEM  FCPP,  the  cell  must  be  supplied  by  a  continuous 
flow  of  hydrogen  to  provide  continuous  energy  output.  Besides 
mechanical  wear  outs  and  control  malfunction,  degeneration  of 
pure  hydrogen  may  lead  to  possible  reduction  in  hydrogen  sup¬ 
ply  to  lead  to  derated  states  of  FCPP  unit.  Since  Q  =  IF  x  wu2 , 
where  Q  is  the  charge  and  wn2  is  the  amount  of  hydrogen 
in  mol;  Q  =  Idc  x  U  then  I^c  =  2 F  x  ^y2-.  Because  qn2  =  ^y2-, 
then  Idc  =  2F  x  qn2,  where  I^c  is  the  total  stack  current  and 
gn2is  the  molar  hydrogen  usage  per  second.  If  the  stack  current 
is  written  in  terms  of  total  stack  power,  then  I^c  =  P/Vs tack,  which 
can  be  rewritten  as 

P  =  2F  -qU2  •  Vstack  (9) 


3.4.  Air  supply  system 


In  a  FCPP,  oxygen  derived  from  air  is  used  to  complete  the 
chemical  reaction.  The  molar  oxygen  usage  per  second,  qo2  of 
the  FCPP  for  a  given  output  power  can  be  computed  as 


P 

Cl02  —  T 77Tt 

4r  Vstack 


(ID 


Eq.  (11)  can  be  adopted  easily  for  air  usage  by  multiplying 
go2with  0.21  since  21%  oxygen  is  available  in  the  air  mix¬ 
ture.  However,  it  is  not  practical  to  consume  all  oxygen  entering 
the  FC  stack.  Hence,  in  real  applications,  the  FCPPs  are  typ¬ 
ically  supplied  with  twice  the  amount  of  air  needed  [8].  For 
this  reason,  the  blower  and  the  compressor  are  sized  so  that 
the  excess  oxygen  ratio  y  is  equal  to  2,  where  y  =  (supplied  rate 
02)/(reacted  rate  of  O2).  Since  more  than  enough  oxygen  is  sup¬ 
plied,  FCPP  output  will  not  be  affected  unless  partial  pressure 
of  oxygen  drops  below  a  certain  threshold  value.  If  oxygen  par¬ 
tial  pressure  falls  below  the  critical  level,  FCPP  will  experience 
oxygen  starvation,  which  may  lead  to  catastrophic  membrane 
failure.  Before  the  system  reaches  this  stage,  an  emergency 
stop  function  will  either  isolate  the  system  from  load  or  shut 
down  the  system.  Hence,  none  of  the  failure  mode  in  airflow 
circulation  leads  to  any  derated  state  but  may  result  in  FCPP 
failure. 


3.5.  Energy  storage  system 

The  flow  rate  of  hydrogen  is  controlled  continuously  in  order 
to  follow  the  electrical  load  variations.  However  depending  upon 
FCPP  type,  flow  rate  adjustment  can  be  achieved  with  a  time 
delay,  ranging  from  a  few  microseconds  to  30  s.  For  this  reason, 
some  type  of  energy  storage  is  needed  for  fast  transient  response 
and  meeting  the  peak  load  requirement.  Among  the  various  types 
of  energy  storage  devices,  batteries  are  the  common  choice, 
but  they  have  relatively  shorter  lifetime  and  need  maintenance. 
Assume  that  the  FCPP  designed  power  is  P  and  system  peak 
load  is  Ppeak-  Thus,  based  on  the  80%  battery  usage,  the  energy 
storage  power  is  Tutorage  =  1 .25  x  CPpeak  —  P)-  The  energy  storage 
capacity  is  measured  in  watt-hours,  expressed  as  W=  Tutorage  x  U 
where  t  represents  the  number  of  hours  used  for  power,  Pstorage- 

If  the  total  power  supply  capacity  of  the  system  falls  below  the 
load  demand  due  to  energy  storage  degradation  and/or  failures, 
it  is  assumed  that  the  excess  load  can  no  longer  be  supplied  and 
ought  to  be  disconnected  from  the  system  by  the  SeMAC  system. 
Hence,  the  percent  reduction  in  the  power  supply  capacity  of  the 
system  due  to  energy  storage  failure  may  be  calculated  as 


From  Eq.  (9),  it  is  evident  that  if  the  maximum  designed  level  of 
power  reduces  by  Agn2,  the  FCPP  output  power  will  decrease 
by  AP  =  2F  •  Agn2  •  Vstack.  As  a  result,  percent  reduction  in 
power  supply  capacity  of  FCPP  due  to  insufficient  hydrogen 


vn7  I(F  +  ^storage)  Ffiemandl 
lS[/c\  =  - 7 - 

DID 

1  1  •'storage 

(P  T  Fstorage)  Fhemancj 


X  100, 


(12) 
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Table  1 

Summary  of  state-space  model  of  FCPP  associated  with  each  system  auxiliaries 


FC  sub-systems 


Failure  severity  FC  system  State-space  representation 


Possible  inadequacy  effects  on  system 
output 


Cooling  system  Inadequacy  down 


Humidification  system  Inadequacy  sown 


Fueling  system  Inadequacy  down 


Air  supply  system  Inadequacy  down 


Energy  storage  system  Inadequacy  down 


Derated  Down 


Derated  derated 


Derated  down 


Down  down 


Derated  derated 


Output  power  will  be  derated  by  RPq% 
due  to  increase  in  temperature  and  inter¬ 
nal  resistance 


Output  power  will  be  derated  by  RPh% 
due  to  increase  in  internal  resistance 
and  decrease  in  partial  pressure  of  the 
reactants 


Output  power  will  be  derated  by  RP^% 
due  to  decrease  in  hydrogen  supply 
and/or  degeneration  of  pure  hydrogen 


Relative  inadequacy  of  air  supply  will 
not  lead  reduction  in  output  power,  but 
may  result  in  complete  system  failure 


Output  power  will  be  derated  by  RPes% 
due  to  FC  stack  and  battery  degradation 


where,  Storage  denotes  the  new  battery  power  due  to  loss  of 
energy  storage  performance.  As  a  consequence,  any  possible 
energy  storage  failure  results  in  only  derated  state.  Other  sys¬ 
tem  components  such  as  stack,  reformer,  power  conditioner,  and 
transformer  can  be  considered  as  essential  components  for  the 
FCPP  power  supply  and  any  of  these  component  failures  may 
bring  the  system  down.  However,  the  stack  has  gradual  per¬ 
formance  deterioration  rather  than  catastrophic  failures.  As  the 
electrodes  and  electrolyte  become  older,  FCPP  output  power 
drops  steadily  with  time.  This  is  more  important  generally  for 
stand-alone  FCPPs.  Table  1  summarizes  the  effects  of  inadequa¬ 
cies/failures  of  FC  sub-systems  on  the  system  output  power  with 
their  state- space  models. 


(14) 


where  T is  the  component  age  (or  the  time  from  last  overhaul),  ft 
the  shape  factor  which  determines  how  the  failure  rate  changes 
with  equipment  age  and  a  is  the  characteristic  time  interval.  If 
P  <  1  the  failure  rates  decreases  with  age,  if  ft  =  1  the  failure  rate 
is  independent  of  age,  and  if  ft  >  1  the  failure  rate  increases  with 
age.  Since  FCPP  technology  is  new  and  there  is  no  publicly  avail¬ 
able  data  associated  with  the  failures  of  FCPP  system,  a  fuzzy 
logic  rule  based  system  is  used  based  on  the  expert  knowledge 
to  determine  the  performance  loss  of  system  auxiliary  compo¬ 
nents.  The  rule  format  for  non-healthy  state  level  is  shown  as 
below: 


3.6.  Estimation  of  transition  probabilities  and  aging  effects  If  input  is  <  component  age  and  maintenance  cycle  and ...  > 


Assume  that  system  auxiliary  component  failure  or  perfor¬ 
mance  loss  takes  place  after  operating  the  FCPP  for  a  while, 
typically  for  5000  h.  Failure  rate  of  the  components  increase 
with  time  until  maintenance  and  drops  to  its  original  value  after 
maintenance.  It  is  assumed  that  maintenance  is  performed  at 
regular  intervals,  and  the  repair  time  and  duration  of  mainte¬ 
nance  increase  when  the  components  become  older.  In  addition, 
the  state  of  health  of  the  auxiliary  components  becomes  worse 
in  proportion  to  component  aging.  The  relationship  described 
above  is  estimated  by  a  Weibull  model,  which  is  a  practical  tool 
for  modeling  component  aging  [9].  The  Weibull  distribution  and 
its  failure  rate  are  defined  as 


f(T)  = 


•  e 


(T/af 


(13) 


then,  output  is  <  degree  of  failure  severity  > 


3. 7.  Reliability  and  EENS  calculation 

In  the  Markov  model,  the  reliability  of  a  system  corresponds 
to  the  summation  of  operating  state  probabilities,  expressed  as 

n 

r  =  (is) 

i=i 

where  Pz  is  a  row  vector  and  represents  the  operating  state  proba¬ 
bilities.  If  the  state  probabilities  corresponding  to  up  and  derated 
states  are  grouped  in  a  vector,  then  Pz  can  be  expressed  as 

P.  _  rpw  PM  Pw  IP^  T>d  p<A 

I  —  L*1  *2  "  •  Vm  lFm+l  Fm+ 2  '  '  '  ^nhxn 


a 


(16) 
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where  represents  the  up  states  and  “d”  represents  the  derated 
states.  Derated  state  implies  that  the  unit  may  not  operate  at  full 
rated  capacity.  Thus  derated  state  probabilities  must  be  reduced 
by  a  factor  i.e.,  correction  vector,  Q,  appropriate  to  the  derated 
state,  given  by 


PEM 

)  (  3- Wire 

LOAD 

(RESIDENTIAL 

HOUSE) 

FUEL  CELL 

CB 

nr 

10  kVA 
120V /240V 

Fig.  7.  Considered  PEM  fuel  cell  system  used  for  reliability  evaluation. 


state  i  can  be  calculated  as 


n EENS);  = 


ft)  (p(0  -  P)dt  +  ff  (/7(f)  -  P)  df 


•8760 


fo  P(0dt 

APi  •  (t i  -  to)  +  AP2  ■  (>3  -  t2) 
*8760 


Jo^P(t)dt 

If  Eq.  (21)  is  generalized,  P(EENS)/  can  be  rewritten  as 

zU  ft(p(t)-p)dt 


(21) 


P(  BENS),  = 


•8760 


Jo  PiO  dt 


(22) 


Ci  =  [l  1  •••  1|C1C2  •••  C„][x„  (17) 

where,  Q  =  l—  RP/%  and  RP/%  represents  the  percentage 
power  reduction  corresponding  to  each  derated  state.  Hence, 
reliability  of  the  system  can  be  calculated  as 

R  =  [P i]\xn  *  [Q'lrcxl  (18) 

As  stated  earlier,  FCPPs  are  subject  to  a  number  of  possible 
outage  and  derated  states  due  to  partial  or  full  failure  of  auxiliary 
components.  Hence,  it  is  essential  to  obtain  an  estimate  of  the 
expected  energy  not  supplied  (EENS)  index,  which  corresponds 
to  the  expected  number  of  watt-hours  curtailed  per  year  due  to 
inadequate  generation,  given  by 

EENS  =  EPNS  x  8760  (19) 

where  EPNS  represents  the  expected  power  not  supplied  per 
hour  at  the  bus  of  interest,  which  is  expressed  as 

V"  (PNS); 

EPNS  =  - -  (20) 

n 

where  (PNS)j  is  the  amount  of  un-served  power  demand  at  the 
FCPP  load  point  associated  with  the  system  state  and  n  is  the 
total  number  of  system  states.  Fig.  6  shows  the  FCPP  behavior 
during  inadequate  power  supply.  From  Fig.  6,  it  is  evident  that 
the  probability  of  EENS,  P(EENS),  at  the  FCPP  load  bus  for 


where,  7  =  1,2,  . . .,  N  is  the  number  of  time  span  (tj+ 1  —  tj)  per 
year  that  the  system  must  curtail  load  due  to  inadequate  FCPP 
generation.  From  the  Markov  model,  the  total  duration  of  derated 
state  that  the  FCPP  may  experience  may  be  obtained  (in  hours)  as 
ti  =  (1  —  Ai)  x  8760,  where  A*  is  the  availability  of  derated  FCPP 
states  that  result  in  lack  of  power  supply. 

4.  Example 

A  methodology  for  reliability  analysis  of  stand-alone  PEM 
FCPPs  is  proposed  and  the  structure  of  the  new  algorithm  has 
been  described  in  earlier  sections.  In  this  section,  the  proposed 
approach  is  applied  for  a  5  kW  PEM  FCPP  that  supplies  a  typical 
residential  house  as  shown  in  Fig.  7. 

The  parameters  of  the  5kW  PEM  fuel  cell  are  shown  in 
Table  2,  where  a  0.0001%  value  is  assumed  for  the  percent  stack 
power  deterioration  per  hour. 

The  residential  power  demand  was  recorded  over  a  800-min 
period  at  a  sampling  rate  of  15  s.  Since  the  capacity  of  the  con¬ 
sidered  PEM  fuel  cell  is  5  kW,  the  balance  must  be  supplied 
by  using  a  standby  energy  storage  system.  The  power  demand 
and  the  power  supplied  to  a  typical  home  through  the  FCPP  and 
a  standby  battery  are  shown  in  Fig.  8(a-c),  respectively,  for  a 
stand-alone  PEM  FCPP. 

From  Fig.  8(c),  it  is  evident  that  the  battery  energy  is  about 
3  kWh.  The  battery  is  recharged  by  the  fuel  cell  when  the  power 
demand  falls  below  the  5  kW  level.  It  is  assumed  that  the  battery 


Table  2 

Five  kilowatt  PEM  fuel  cell  parameters 


Parameter 

Value 

Parameter 

Value 

Cell  number  (No) 

Active  area  of  a  cell 

Universal  gas  constant  ( R ) 

Temperature  of  the  fuel  cell  stack  ( T) 
Faraday’s  constant  (F) 

Hydrogen  partial  pressure  (puA 

Water  partial  pressure  (pibo) 

100 

100  cm2 

8. 1345  J  moU1  K-1 

353  K 

96485  C  moU1 

1.087  atm 

0.464  atm 

Oxygen  partial  pressure  (po2) 

Charge  transfer  coefficient  of  the  electrodes  (a) 

Current  of  the  FC  stack  (7dc) 

Limiting  current  of  FC  stack  (Ium) 

Exchange  current  density  (7o) 

Empirical  coefficient  for  concentration  voltage  (c) 

Stack  internal  resistance  (Rfc) 

2.17  atm 

0.5 

94.69  A 

105  A 

10— 6  912  A  cm-2 

0.0147 

0.00303  Q 
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Fig.  8.  (a)  Typical  residential  power  demand,  (b)  fuel  cell  power  and  (c)  battery  power  to  the  load  vs.  time. 


is  replaced  periodically  (5  years  interval),  which  is  defined  as 
the  point  where  the  capacity  has  declined  to  80%  of  the  nominal 
value. 

The  proposed  FCPP  modeling  incorporates  the  component 
age  and  maintenance  cycle  (input  values)  and  determines  the 
degree  of  failure  and  failure  rates  via  membership  functions  and 
Weibull  distribution  function.  It  is  assumed  that  the  possible 
wear  out  period  starts  after  5000  h  of  operation  and  continues  to 
maintenance  time,  typically  1  year.  This  process  is  repeated  at 
regular  maintenance  intervals  and  the  repair  times  and  the  dura¬ 
tion  of  maintenance  increase  at  each  cycle  as  the  state  of  health 
of  the  auxiliary  components  decreases.  The  states  of  health  of 
subsystems  are  represented  individually  by  the  triangular  mem¬ 
bership  functions  as  shown  in  Fig.  9,  which  are  developed  by 
using  expert  knowledge.  In  Fig.  9,  since  maintenance  is  supposed 
to  be  done  at  regular  intervals,  component  age  gives  direct  infor¬ 
mation  about  the  number  of  maintenance  that  has  been  done  so 
far.  Therefore,  input  membership  function  is  chosen  for  only 
component  aging. 

Rules  used  in  the  fuzzy  inference  system  are  listed  below: 

•  If  age  is  <A\  (very  young)>,  then  reduction  is  <B\  (very 
healthy)> 

•  If  age  is  <A2  (young)>,  then  reduction  is  <B 2  (healthy )> 

•  If  age  is  <A3  (medium)>,  then  Reduction  is  <B 3  (medium)> 

•  If  age  is  <A4  (old)>,  then  reduction  is  <B 4  (worn)> 

•  If  age  is  <A5  (very  old)>,  then  reduction  is  <B$  (very  wom)> 


0  4  8  12  16  20  24 

(a)  Age  [year] 


These  rules  show  the  functional  relationship  between  component 
age  and  reductions  in  auxiliary  performance.  Fig.  10  depicts  the 
state-space  modeling  of  a  5  kW  PEM  FCPP  generation  for  a  5- 
year  operational  period.  It  is  assumed  that  the  stack  and  battery 
performance  degrades  gradually  with  time.  Hence,  the  stack  and 
battery  energy  generation  sets  are  merged.  The  transition  prob¬ 
abilities  in  the  model  are  estimated  for  each  subsystem  based  on 
the  Weibull  distribution. 

5.  Simulation  results 

The  simulation  of  FCPP  reliability  variations  using  the  pro¬ 
posed  approach  is  carried  out  for  different  operational  lifetimes 
as  shown  in  Fig.  11.  The  instantaneous  FCPP  reliability  varia¬ 
tions  for  operational  period  of  1-5,  and  6-10  years  are  shown 
in  Fig.  11(a)  and  (b),  respectively.  Fig.  11  shows  the  FCPP  reli¬ 
ability  as  a  function  of  stack  and  battery  degradation  indicating 
that  the  reliability  decreases  with  time. 

Fig.  12  shows  the  steady-state  reliability  as  well  as  cumulative 
EENS  variation  of  FCPP  as  a  function  of  time.  In  Fig.  12,  the 
solid  line  represents  the  reliability  of  FCPP,  while  the  dashed 
line  represents  the  cumulative  EENS  values,  which  are  given 
for  the  first  and  second  5 -year  period  of  operation.  The  results 
of  Fig.  12(b)  is  obtained  assuming  that  the  battery  is  replaced  at 
the  end  of  5 -year  operational  period. 

It  is  evident  from  Fig.  12  that  the  reliability  level  drops  faster 
with  time.  The  reliability  increased  considerably  in  Fig.  12(b) 


5  10  15  20  25 


(b)  Auxiliary  State  of  Health  Reduction  [%] 


Fig.  9.  An  example  of  input  (a)  and  output  membership  (b)  functions  for  state  of  health  of  fuel  cell  auxiliaries. 
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Fig.  10.  The  state-space  model  ofa  5  kWPEMFCPP  generating  unit  for  the  fifth  year  operation,  where  Ai  =  1.1423  x  10  4,  A.  2  =  1.1422  x  10  4,  A.3  =  1 . 1421  x  10  4, 
A. 4  =  1.1419  x  10~4  [failure/hour]  and  /x\  =0.3088, 1x2  =  0.3388,  /x 3  =  0.3758,  p4  =  0.2857  [repair/hour]. 


Fig.  11.  Instantaneous  FC  system  reliability  (a)  from  1  to  5  years,  and  (b)  from  6  to  10  years,  respectively,  from  upper  to  lower  lines. 


Fig.  12.  Steady-state  reliability  and  cumulative  EENS  variation  of  FCPP  as  a  function  of  time  for  the:  (a)  first  5-year  lifetime,  and  (b)  second  5-year  lifetime  with 
battery  replacement  (solid  line:  reliability  variation,  dashed  line:  EENS  variation). 
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Fig.  13.  Steady-state  reliability  and  cumulative  EENS  variation  of  FC  with  load 
management  system  for  10-year  operational  period  (solid  line:  reliability  varia¬ 
tion,  dashed  line:  EENS  variation). 

due  to  battery  replacement.  However,  2  years  after  battery 
renewal,  a  noticeable  reliability  decrease  started  due  to  stack 
degradation  from  the  first  5 -year  operation  period. 

In  general,  a  FCPP  may  be  unavailable  for  a  significant  dura¬ 
tion  due  to  lack  of  generation  capacity  resulting  from  stack 
and  battery  degradation.  However,  this  negative  impact  can 
be  minimized  by  utiliziing  a  smart  energy  management  con¬ 
trol  (SEMaC)  system  in  stand-alone  FCPPs.  A  SEMaC  system 
can  intelligently  schedule  the  power  demand  based  on  the  user 
load  profile  and  priority  of  needs,  especially  to  tackle  excess 
energy  requirements.  The  steady-state  reliability  and  cumula¬ 
tive  EENS  variation  of  the  FCPP  as  a  function  of  operational 
age  is  depicted  in  Fig.  13.  It  is  evident  from  Figs.  12  and  13  that 
the  FCPP  availability  can  be  increased  significantly  by  imple¬ 
menting  the  SEMaC  system.  For  example,  using  SEMaC  based 
load  management  system,  a  customer  may  experience  as  low  as 
90  min  unavailability  (instead  of  lOh)  during  the  first  year  and 
4  h  unavailability  (instead  of  5  days)  in  the  5th  year  of  FCPP 
operation. 

6.  Conclusions 

In  this  paper,  the  operational  benefits  from  load  management 
are  evaluated  for  a  stand-alone  PEM  FCPP  by  implementing 
a  novel  FCPP  reliability  assessment  technique.  The  proposed 


technique  has  been  found  to  be  an  effective  tool  for  model¬ 
ing  the  reliability  of  FCPP  and  associated  systems.  The  results 
obtained  using  the  proposed  technique  confirms  that  the  reliabil¬ 
ity  of  FCPPs  drops  steadily  with  time  as  the  components  become 
older.  However,  FCPP  reliability  can  be  considerably  improved 
by  utilizing  the  SEMaC  based  load  management  and  control  sys¬ 
tem,  which  also  increases  the  competitiveness  and  market  value 
of  the  stand-alone  FCPPs.  The  proposed  FCPP  reliability  eval¬ 
uation  method  can  be  used  as  a  tool  to  facilitate  maintenance 
scheduling  of  FCPPs  such  as  battery  replacement.  Apart  from 
these,  bipolar  plates  and  membrane  electrode  assemblies  can  be 
renewed  in  order  to  lengthen  the  lifetime  of  FCPPs.  The  sched¬ 
ule  for  this  replacement  can  be  effectively  decided  by  utilizing 
the  proposed  approach. 
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